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Abstract 

This paper is concerned with the restriction of interval exchange transformations 
to algebraic number fields, which leads to maps on lattices. We characterize renormal- 
izability arithmetically, and study its relationships with a geometrical quantity that 
we call the drift vector. We exhibit some examples of renormalizable IET with zero 
and non-zero drift vector and carry out some investigations of their properties. In 
particular we look for evidence of the finite decomposition property on a family of 
IETs extending the example studied in [18]. 
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1 Introduction 



An interval exchange transformations (IET) is a one-dimensional map E of an interval, 
consisting of a rearrangement of a partition of that interval. Let {Jlj} be an iV-elements 
partition of £1 = [0, 1), and let ir be a permutation of the integers 1,... ,N. We define 
E as the map that rearranges the intervals according to the permutation tt, preserving 
orientation. This procedure defines N translation vectors Tj 

E(x) = x + n x^VLi. (1) 

We assume that this permutation is irreducible, that is, ir{l, . . . , k} = {l,...,k} if and 
only if k = N. 

There is substantial literature on the subject. An early theorem of Keane [13] states 
that if the permutation is irreducible and the set of lengths rationally independent, then 
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two discontinuity points never lie on the same orbit and E is minimal. A stronger results 
was later obtained by Masur and by Veech [20,25]: given an irreducible permutation, for 
almost all set of lengths, the associated interval exchange map is minimal and uniquely 
ergodic. 

Boshernitzan [8] gave the following ergodicity criterion: if we let e& be the shortest 
length of an interval of the partition of E k , then if 

lim ket > 

E is uniquely ergodic. A result of Katok [12] ensures that interval exchanges cannot be 
mixing. However, a recent result by Avila and Forni [3] shows that, for any permutation 
which is not a rotation (meaning that ^ i + 1 (mod N)), the interval exchange map 
associated to Lebesgue-almost all set of lengths is weakly-mixing. 

Here we consider a uniquely ergodic interval exchange map E, and we assume that 
the lengths of the sub-intervals fli generate an algebraic extension K of degree n of the 
rationals. The module of E is the set of all linear integral combinations of the sub-interval 
lengths, namely 

N 

M = {^mil^l : nii £ z} (2) 

i=l 

and the rank of E is the rank of its module. Algebraic IETs constitute a countable, hence 
exceptional, family. Little is known about these maps, the only substantial result being 
the Boshernitzan and Carroll theorem [9]: if an IET is defined over a quadratic field, then, 
up to rescaling, the number of induced maps on sub-intervals is finite. In short, quadratic 
IETs can be renormalized. In general, interval exchanges over algebraic number fields of 
higher degree do not have this nice property. If they do, then their rank must be greater 
than two. 

The code u> = uj{x) of a point x G fl is constructed with respect to the natural partition 

ujj = i 44> E j ~ 1 x € 

Let ei, . . . , ejv be the canonical basis of Z N . We associate to each word uj its sub-word 
of length k, and the sequence s(lo^) in 7L N defined as follows 

k 

S (u,«) = J>^ A; = 1,2,.... (3) 

j=i 

We call it the staircase function. If we denote by the number of occurrences of i in 

then we have 

N 

s(^) = J2\^%e t - 
i=i 

This expression may be rewritten as follows 

s(J k) ) = kA + D k . 
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Here A is the vector of the lengths of the intervals 

A=(|0 1 |,...,|O iV |) (4) 

and D k is the vector of the local discrepancies of the sequence u in each cylinder (i) or, 
equivalently, of the sequence (Ex) in each intervals £2j 

D k = (D k {x,n 1 ),...,D k (x,n N )) (5) 

and 

k 

D k (x,n i ) = J2^ i (E j=l x)-k\n i \ = \J%-k\ni\, 

3=1 

where 1^ is the characteristic function of the interval Oj. By Birkhoff's ergodic theorem 
and using minimality we have ||.Dfc|| = o(k); the staircase of each point thus has the 
principal direction A with sub-linear perturbations. 

The theory of the deviation from the ergodic averages in IET is closely related to 
the theory of geodesic flows on flat-surfaces, for which strong results have been given by 
Zorich [27,28]. Specifically, we suppose that the IET can be described as a section of a 
geodesic flow on a flat surface of genus g. Then, for almost all such IETs, there is a nested 
sequence of subspaces 

Hi c Ti 2 c • • • c n g c R N 

with Tii being the one-dimensional subspace spanned by A, and each Tii being z-dimensional, 
such that, for any / £ Tif- with ||/|| = 1 

r log(g(a/ fc >),/) 

llHl SUp : = V{ . 

fc^oo log k 

(The symbol ( , ) denotes the standard inner product.) This result says that the deviations 
of the staircase function from the ergodic means follow some power laws, the exponents 
of which are given by the Lyapunov exponents z^, depending only on the so-called Rauzy 
class the IET belongs to. 

The aim of this paper to introduce a lattice representation of algebraic IETs, following 
a method developed for higher dimensional piecewise isometries [14,16], and analogous to 
the representation used in [18]. Since the translations Tj are expressed as integral linear 
combinations of the lengths, for any point £ in K D [0, 1) the set £ + M is invariant under 
E and contains the orbit of £. We then identify £ + M with the lattice = Z m , where 
m < n is the rank of M. This construction may be seen as a projection of the lattice Z N 
in which the staircase is built onto the lattice L where the dynamics takes place. In this 
way we obtain a discrete geometric representation of the IET, restricted to a dense set of 
algebraic points. 

The content of this paper is the following: 

In section[2]we introduce a module-theoretic description of an algebraic IET and define 
a fundamental dynamical invariant, the drift vector S. We then show that the drift vector 
is a representation of the classical Sah-Arnoux-Fathi invariant (proposition [5]). Most of 
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the constructs developed in this section refer to the case of maximal rank; this restriction 
is justified by the fact that renormalizable IETs have this property (see theorem [1] below). 

In section [3] we deal with renormalizability of algebraic IETs and we describe the role 
of the drift vector in the renormalization process. We then prove the main result of this 
section, which appears as theorem [6] and proposition [HJ (For background on substitution 
dynamical systems see [22].) 

Theorem 1 Let E be a uniquely ergodic interval exchange transformation on N intervals, 
which is renormalizable with scaling ratio p. Then p is a unit in the ring Z[/3], where 
(3 = p~ l is the largest eigenvalue of the incidence matrix of the primitive substitution 
associated with E. The map E is defined over Q(p), its rank n is maximal and equal to 
the degree of p, and 

(i) if the drift vector is non-zero, then n is even and n < N; 
(ii) if the drift vector is zero, then 2n < N. 

Subsequently, we introduce the recursive tiling property and formalism (Vershik cod- 
ing), and give estimates for the asymptotic behaviour of points with eventually periodic 
code in the scaling dynamics (propositions [9] and [TT1) . We also give a condition for these 
points to be of positive density on the lattice (proposition [TO]) . 

Section [4] is devoted to the analysis of specific examples. Using Rauzy induction, 
we construct a renormalizable IET of degree four over four intervals, with non-zero drift 
vector. We then construct an infinite family of cubic renormalizable IET with zero drift- 
vector, which belong to the so-called Arnoux- Rauzy family [6]. We also perform numerical 
investigations which give strong evidences that some of these examples share the finite 
decomposition property encountered in [18]. This property, which we establish in one case 
(proposition [T3|) is described as follows. The field K may be represented as the disjoint 
union of invariant layers £ + M, for a suitable set of representative points £ € K. An 
algebraic IET has the finite decomposition property if the restriction of each layer to the 
unit interval is the union of finitely many orbits. 

These results show connections between the finite decomposition property of a renor- 
malizable IET and the Pisot property of the (reciprocal of the) scaling constant. This 
problem deserves further investigation. 

Acknowledgements. This research was supported by EPSRC grant No GR/S62802/01. 
One of the authors (GP) wants to thank S. Ferenczi, X. Bressaud and P. Arnoux for some 
enlightening discussions. 
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2 Lattice dynamics and drift vector 



In this section we consider algebraic IETs with maximal rank, for which an interesting 
arithmetical theory arises. In the next section we shall see that renormalizable IETs have 
this property. 

2.1 Lattices 

We choose a basis for M. as a Z-module, and consider the translation vectors T{ £ M. 
defined by E — cf. equation ([I]). Let be the map from M. in L sending the basis 
elements of A4 onto the canonical basis of L. We have two lattices, one is the symbolic 
lattice Z N in which the staircase was defined; the other is L, it has dimension n and 
represents the module in which the dynamics takes place. We call n the projection from 
the former in the latter; this is a linear map such that 

vr(ei) = Vi Vi := 0(rj) (6) 

which is represented by an n x N matrix. We denote by tp = (ftoEocj) -1 the map conjugate 
to E on the lattice L. 

As we have a map of the unit interval, the dynamics of ip does not take place on all of 
L, but only on the invariant slab of points z £ L such that 6 [0, 1). We are then 

led to consider the set M n [0, 1) for which we require a suitable representation. Let 

O = {C G K : CM C M} (7) 

where K is the smallest field containing A4. Then O is an ordei0 in K. Let v\, . . . , v n 
be a Z-basis for M, and let d be the smallest positive integer such that dv\~ G O, for 
k = 1, . . . ,re. Then J = dM. is an ideal of O [10, section 4.7]; the normalized generators 
di^k, k = 1, . . . ,n, are algebraic integers which form a Z-basis for J. By construction, Ai 
is also a module over 0; furthermore O is the largest subring of K with this property. 
Now 5^i=i = 1 € an d so Z C M. We can form the quotient 

M! :=M/Z 

which is a Z-module, which we shall represent asMn [0, 1). Let j be the smallest positive 
integer in J: then the rational numbers in M are the multiples of j/d, and hence 

(QnM)/z = i-z/z^z/bz b = — -f— , (8) 

d gcd(ii,j) 

that is, there are b congruence classes of rational numbers in A4', while the non-rational 
part is freely generated. Thus we have a Z-module isomorphism M' — * L' where 

L' = (Z/6Z) e Z n_1 (9) 
A finitely-generated subring of K, which contains 1 and has maximal rank. 
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which characterizes the additive structure associated with the map E. Accordingly, we 
choose for Q G Ai' the following representation 

1 n 

C = 2 mk u 'k m k eZ v' k G J, v[= j, (10) 

k=l 

where the vL are obtained from the normalized generators dv k via a unimodular transfor- 
mation. The isomorphism ([9]) then becomes 

4> : £ i-> (m (mod 6), mi, . . . ,m n _i). (11) 



To represent the dynamics in the field K — the extension of the rational generated 
by the intervals' lengths — we use the Z-basis {v k ) for Ai, and any point x in K can be 
written as 

n 

x = ^r fc ^ fc for r k G Q. 
k=l 

By definition, E preserves Ai and for any x G K its orbit is confined in x + M. We have 
a "layered" representation of K and we regard it as the union of disjoint invariant sets of 
the form 

K=[j{i + M) (12) 
for a chosen set of representative H. Following [18], we choose 

n 

Z = {j2 r kVk, r fc GQn[0,l)}. (13) 
fc=i 

The dynamics of E acts on each layer separately, which in each case may be represented 
as the lattice L, with a ^-dependent dynamics on it. 

In general, the closer two points get in A4, the further away their images in the lattice 
L will be. Indeed, let x G A4 correspond to the lattice point z = (z%, . . . , z n ) 

n 

x = J2z k u k e M. (14) 

k=l 

If the coefficients z k are bounded, then there are only finitely many possibilities for x, and 
hence x is bounded away from zero. Thus if \x\ — > then ||z|| — ► oo, and one is interested 
in the relation between these limits. 

We recall that height of a polynomial is the maximum absolute value of the coefficients, 
and the height of an algebraic number is the height of its minimal polynomial. In the 
particular case when A4 takes the form Ai = Z[A] for some algebraic number A of degree 
n, El Liouville's inequality gives the following bound. 

2 Z[A] is the Z-module generated by 1, A, A 2 , ... , A' 1-1 
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Proposition 2 If A4 = Z[A] for some algebraic number A of degree n, then for every 
x £ M we have 

\x\ \\z\\ c > e^-V z = <f>(x) c = degA + logtf(A) (15) 

where H(X) is the height of X. 

Proof: We have (cf. [24]) that for an integer polynomial / with height < H and degree 
d, for all algebraic number A we have either 

|/(A)| > e -< d +^ H ) 

or A is a root of /. The constant is c = degA + H(X). In our context, the polynomial / 
has coefficients Zi, and A is the generator of Ai, that is, the basis vectors = \ k ~ 1 . Then 
H = \ \z\\ for the max-norm and 

|/(A)| = |x|>e- c ^- 1 )||z||- c , 

as desired. □ 



2.2 Drift vector 

Let Vi be the image on the lattice of the translation vector n — see equation (|6|). The drift 
vector S, defined as 

N 

S = \Oi\vi = tt(A) (16) 

i=l 

gives the dominant direction of the dynamics over the lattice. The drift vector is related 
to a known invariant of IET theory, the so-called Sah-Arnoux-Fathi invariant, as we shall 
see in section 12.41 

Let the point x G M have itinerary u>, and let z = 4>{x) be the corresponding lattice 
point. Then the orbit of z may be expressed in terms of the drift vector as follows 

z + vrs(w (fc) ) = z + kS + irD k = z + k(S + a(k)) (17) 

where the staircase function s and the discrepancies were defined in equations (|3]) and 
([5]), respectively, and a(k) = nD^/k. Because \Df.\ = o(k), the orbits of any lattice map 
with nonzero drift vector will be dominated by a linear drift in the direction of S. 



2.3 Density 

The dynamics in L takes place on a subset of zero density. We can nonetheless have a 
meaningful measure of the "size" of the orbits by defining a modified density function on 
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the reduced lattice I/, which results from restriction to the unit interval. With reference 
to dHj) and ([!]), we represent a subset A of M! as 

A = Ao + '-' + Aj,-! AinJ2 A k = {n} A t nq= l -{modZ). 



b 



Letting Aj = (j)(Aj), we define 



D ( Aj ^ = k^L b(2k) n ~ 1 *- m ° d ^' mi ' ' ' ' G A?> l mi l ^ k ^¥= 0}- 

Note that mo (mod 6) is fixed, by construction. If the limit above exists for all j, we define 
the density of A as 

6-1 

Z>(A) = J>(Ai) 

i=o 

which is a (non countably additive) probability measure on L'. 

The density T> is invariant under an unimodular transformation U. The density T> 
measures the proportion of the points lying inside an hyper-cube Ct centred at the origin, 
with side 2k, as k becomes large. Considering D(UA) is equivalent to consider D'(A) where 
V measures the limiting proportion of the points lying inside the rhombus U~ l Ck, also 
centred on the origin. As for all k there exists k! and k" such that U~ 1 Ck> CC^C U~ 1 Ck" 
we deduce that V(A) = V'(A) = V(UA). 



Proposition 3 Let I C [0, 1) be a finite union of intervals, and let A = <p(M n /). Then 
V( A) exists and is equal to the Lebesgue measure of I. 



Proof. We begin with the case in which / is a single interval. Let m = (mo, • • • , rn n ^\) 
let x m = _1 (m), and let 

C k = —— - — rt){m G A : I mil < k\ k = 1,2, 

Now let mQ,ni2, ■ ■ ■ , m n _i be fixed integers, and define 

Cm ,fc,m 2 ,...,m n _i = ( m > m j m 2 , • • • , ^n-l) : \m\ < k, X mo ,m,m 2 ,...,m„_i G 1} 

then by Weyl's criterion 

Cmo,fc,m 2 ,...,m„_i - |I| ->■ when k -»• oo. 

We then have 

j 6-1 fc 

^ k ~ b(2k) n ~ 2 ^mo,k,m2,—,m n -i 
mo=0 m2,...,m„_i=- fe 

and 

. 6-1 fc 

^ fc = b(2k) n - 2 ^ ^ ( e m ,fc,m 2 v,m„_i + 1-^1) • 

mo=0 m2,...,m„_i=- fc 
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where e m0) js !)m2) ... )mn _ 1 = C mo ,fc,m 2 ,...,m„_i - |/|. Finally 

j 6-1 fe 

Cfc = l-^l + U2k) n - 2 ^ ^ e m ,fe,m 2 ,...,m„_i 

mo=0 m2,...,m„_i=- k 

the last term is a (n — 2)-dimensional Cesaro sum and thus converges to zero when k — > do, 
leading to 25(A) = |/|. The argument above is easily extended to the case of a finite union 
of intervals. □ 

Next we have 



Proposition 4 If the drift vector is non-zero and if n > 2, then the density of any orbit 
on the lattice L' is zero. 



Proof. This is an easy consequence of formula (|17p . We have 



\\t(j k (z) -z\\ = \\Sk + ka(k)\\ > k 



\a(k)\\ 



> fcmin 

k 



\a(k)\\ 



This latter minimum is attained because 1 > ||a(/c)|| — > 0. Thus the number of points 
inside the box centred at the origin and with side 2k increases at most linearly with k. So, 
if the degree of the IET is greater than 2, then the lattice L' has rank greater than one, 
leading to zero density. □ 



In particular, an interval exchange transformation can have the finite decomposition 
property — each of its layers is the union of a finite number of orbits — only if it has 
zero drift. The actual value of the drift vector depends on the choice of the basis of A4; 
changing the basis is equivalent to apply a unimodular transformation on S which would 
indeed change its direction and magnitude. Nevertheless, we will see that the relevant 
property of S is mainly its nullity. 



2.4 The Sah-Arnoux-Fathi invariant 



In this section we regard R as a vector space over Q, and consider the tensor product 
R® M of Q- vector spaces. To a general iV-interval exchange E with set of length \Vti\ and 
translation vectors r.;, Arnoux [5] associates an element S of M ® R. The quantity S can 
be computed by the formula 

N 

s = i f2ii ® Tj. 

i=l 

The Sah-Arnoux-Fathi (SAF) invariant is a non-trivial conjugacy invariant; it is not a 
complete invariant but is in some sense the best possible one. The set of all the IETs 
constitute a group under the composition law. The map associating to each IET its 
invariant S is a homomorphism from the group of IETs into R ® R, with the smallest 
possible kernel. For further details and discussions about the algebraic significance of S 
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see [5,26]. From these references we also know that S is an induction invariant for minimal 
transformations, i.e., S is the same computed on any induced maps of a given minimal 
interval exchange transformation. If E is periodic, then its invariant is zero. The converse 
is true for an interval exchange of three sub- intervals, but not for more subintervals, as 
shown in [5]. The Arnoux-Yoccoz maps [4,7] provide an example of this phenomenon: 
they have 7 sub-intervals and zero invariant, but they are not periodic, in fact they are 
minimal and uniquely ergo die. 

In general, R(8>M is infinite dimensional, but in our case the lengths and the translations 
belong to K, which is a vector space over Q of dimension n. In fact we can even restrict 
ourselves to the Z-module M. . The quantity S should then be in M. ®z M. (the subscript 
Z emphasizes that the tensorial product is taken between Z-modules). As above, we fix a 
Z-basis V\ y . . . , v n for M.. We define a linear map V : M. (&i M. — > M. n by specifying its 
action on this basis as follows 

V : Ui <g> Vj i — > (0, . . . , Ui, . . . , 0) . 

v „ ' 

Vi has position j 

The map V is bijective and is an isomorphism of Z-modules. We have 

N 



v(s) = Y / v(\n i \®T i 



1=1 



Next we express the translation vectors ti and lengths with respect to the module 
basis 



n 



fc=i j=i 



so that 



j,k=i 

Finally, recalling that = 0(Tj) are the images of the translation vector on L (equation 
flEj)), we have 

i=l j=l 3=1 

N N 



i=i i=i 



which shows that the drift vector S is a representation of the SAF invariant in our specific 
context. 



The condition S = is quite strong; in fact it implies n rational dependencies between 
the lengths of the intervals. Note that if we have only two intervals, the two lengths are 
of the form Q\ and 1 — fli and the condition 5 = would imply \£li\k + (1 — |fii|)Z = 
for two integers k and /, which is impossible since is irrational. In fact a two-interval 
exchange is a pure rotation of the circle, which gives a pure translation on the lattice. 
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More generally, for N intervals, we have N lengths G M, which we assumed to have 
rank n. For S to be zero, we must have 



for some integers x\. Then a necessary condition for 5 = is n < N — 1. 
The following result summarizes the considerations above. 

Proposition 5 The SAF invariant of an algebraic IET is zero if and only if its drift 
vector is zero. If the degree n and the number of intervals N are such that n > N — 1, 
then the drift vector, and hence the SAF invariant, are non-zero. 

By the discussion above, the Arnoux-Yoccoz [4, 7] and Arnoux-Rauzy [6] families of 
IET have zero SAF invariant and thus zero drift vector. 



3 Renormalizable cases 

3.1 Algebraic consequences of renormalizability 

In the following, we assume E to be renormalizable. An IET of N intervals is renormaliz- 
able if it has the same dynamics at each scale level, meaning that there exists an interval 
I C [0, 1) with the property that the induced map on / is conjugated to the IET itself via 
a simple rescaling p. Thus the induced map on / is an IET of N intervals whose lengths 
are p\£li\ with a common scaling factor p. In this section, we prove the following theorem. 

Theorem 6 Let E be uniquely ergodic interval exchange transformation on N intervals, 
which is renormalizable with scaling ratio p. Then p is a unit of degree < N in the ring 
Z[/3], where [3 = p~ x is the largest eigenvalue of the incidence matrix of the primitive 
substitution associated with E. Furthermore, E is defined over Q(p) and it has maximal 
rank. 



Proof. By construction, the scaled lengths can be expressed as integral linear combina- 
tions of the lengths If we represent these linear combinations as an N by N integral 
matrix B, then we see that, by definition, the lengths vector A is an eigenvector of B with 
eigenvalue p. This tells us that p is an algebraic integer, and if n is its degree, then n < N . 
Since B ■ A = pA and A is a probability vector, then 



N 



N-l 




N 



N 



■ A)i = p 



i=i 



i=i 



ii 



We see that pAi C Ai and p G Ai. Now, Ai contains 1, and repeated scaling show that 
it also contains all positive powers of p. Hence Ai contains the ring Z[jo]. It follows that 
the rank of Ai is at least as large as the degree n of p. 



A renormalizable and uniquely ergodic interval exchange is conjugated to a primitive 
substitution a. Furthermore, the atoms of the induced map on I have their iterates cover 
the whole interval f2, lest there would be either some periodic cells or some different ergodic 
component. We conclude that 

JV 

^p k \a k {i)\\Sli\ = l k>0. 

i=l 

Because the substitution is primitive, for all i the quantity |<J fc (i)| grows like p k , the powers 
of the Perron eigenvalue of the incidence matrix M a of a, and we have 

Because both p and (5 are algebraic integers, they are units. 



Next we show that A is the positive eigenvector of M a associated with (3. Considering 
the fixed point u of a, let j be the unique integer with the property that a k (j) is an 
approximation to lo. We find 



\ \n N \ 

For all i, we have 




lim , 



- k+1 m, 



lim 

k— >oo 



1 



W k (j)\ 



V w k+1 m 



N 



W k (j)\ 



W k (j)\ W k +i(j)\ 



P\Sli\. 



Since M a A = (5 A, for the same reason as above we have (3 G Ai, and as (3AA C Ai and 
fj~ l M = pAi C Ai, we have j3M = pM = Ai. We see that the vector of the lengths is 
the solution of the system of equations 

JV 

M a z = (3z and Zi = 1 z = (z±, . . . , zn). 

i=i 

This solution is unique. Indeed, by the Perron- Frobenius theorem, dimker(M CT — (31) = 1, 
and the solution belongs to Q(/?) = Q(p). Then Ai C Q(p), and therefore the rank of Ai 
is at most n. But then the rank is exactly n, namely, it's maximal. □ 

The proof above shows that Ai is a Z[p]-module, and hence Z[p] is a subring of the 
order O — see equation ([7J and following remarks. 

Let R be the matrix representing the action of multiplication by p on the lattice L. 
Since p is a unit, the matrix R is unimodular; moreover, from linear algebra and the 
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fact that p has degree n, we have that the characteristic polynomial of R is equal to the 
minimal polynomial of p. In particular, p is an eigenvalue of R. 



The codes of the orbits from a scale to the next can be found by mean of the substitution 
a. If a point x belongs to the atom i, we have 

E^(px)=pE(x), (18) 

or on the lattice 

^^(Rz) = Rnj}{z). (19) 
Letting the projection tt be as in (JH), we have 

\cr(i)\ N N 

j=i j=i j=i 

N 

= Rz + nM a \i\j e j = R z + ^MaCi. 
3=1 

By equation (JT9J) , we then have 

Rip(z) = Rz + Rvi = Rz + TrM a ei 

and since irei = Vi 

Rirei = irM a ei i = 1, . . . , N. 

Finally 



Lemma 7 The following diagram commutes 







7L N 


I* 




I- 


L 


R 


L 



As a corollary, we see that if X is an eigenvector of M a with eigenvalue C and X is not in 
the kernel of tt, then irX is an eigenvector of R with the same eigenvalue 

irM a X = (ttX = RttX. 

In addition, the matrix M a preserves ker-zr. Indeed if X £ ker7r, then again R*kX = = 
TxM a X and M a X £ ker-zr. In particular, the characteristic polynomial of M a has a factor 
of degree N — n. 

We already know that A is an eigenvector of M a with the Perron eigenvalue (3. We thus 
have two possibilities: i) ttA ^ (the drift vector is non-zero) and irA is an eigenvector 
of R with eigenvalue 0; ii) A £ ker7r and (5 is not part of the spectrum of R. This is 
because the spectrum of R is contained in the spectrum of M a and the Perron eigenvalue 
is simple. 
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3.2 Connections with drift vector 

When the drift vector S is non-zero, the algebraic number (3 is in the spectrum of R. This 
means that p and (5 = p^ 1 are algebraic conjugates, namely roots of the same irreducible 
polynomial p(x) of degree n. In fact, all eigenvalues of R are algebraic conjugates and all 
are units. 

Let now p(x) = x deep p(l/x) be the reciprocal polynomial of p — the polynomial whose 
roots are the reciprocal of the roots of p. Because p is irreducible, so is p, and since (3 is 
a root of both p and p, it follows that p = p. So the roots of p regroup into reciprocal 
pairs, and since these roots cannot include ±1, the degree n of p must be even. Any 
self-reciprocal polynomial p of even degree takes the form (see for instance [2]) 

p{x) = x n/2 h(x + x- 1 ) 

for some an integral polynomial h of degree n/2. The polynomial h is irreducible because 
p is. So if C is root of p, then £ = £ + £ -1 is root of h, and since £ 2 — ££ + 1 = 0, we see 
that £ lies in a quadratic extension of the algebraic number field generated by £. 

Now p is the characteristic polynomial of both R and R -1 , which have the same 
spectrum. Because p is a minimal polynomial, its roots all have multiplicity one and 
therefore R and R~ l are diagonalizable (see, for instance, [11, §3.10]). It then follows that 
R and R^ 1 are conjugate, indeed integrally conjugate. That is, there exists an integral 
matrix G such that 

R- 1 = GRG- 1 

showing that R has time-reversal symmetry [15]. Then G is an involution, exchanging two 
subspaces of dimension n/2. 

When the drift vector S is zero, the eigenvalue (5 of M a is not in the spectrum of R. 
Then the characteristic polynomial P{x) of M a admits the factorization 

P{x) = p{x)p{x)q{x) (20) 

where the polynomials p and p are irreducible of degree n, while q is a polynomial of degree 
N-2n. 

Summing up the results of this section we have 

Proposition 8 Consider a uniquely ergodic algebraic IET on N intervals, degree n and 
maximal rank. 

i) If the drift vector is non-zero and n is odd, then renormalization is not possible; 

ii) if the drift vector is zero then 2n < N. 
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Note that if we have only two intervals, the two lengths are of the form J7i and 1 — fii 
and the condition 5 = would imply (1 — |^i|)^ = for two integers k and I which 

is impossible as is irrational. In the case when n = N = 2, the drift vector must 
be non-zero and yet the map is known to be renormalizable, from the Boshernitzan and 
Carroll theorem [9]. 

3.3 Recursive tiling property 

We already introduced the standard coding, with respect to the partition of the interval 
n. In [14,16,18] and elsewhere, a different type of coding is also used, which is based on 
the so-called recursive tiling property. As we mentioned earlier (in the proof of theorem 
[6J, for each k, the interval Q splits into YliLx & k (i) sub- intervals, each corresponding to a 
fc-cylinder in the symbolic space. For a given point x, knowing for each k in which cylinder 
x lies is the essence of this new coding, which, for reasons explained later, is also called 
the Vershik coding. 

We consider the substitution a introduced above, using the notation lo 1 = cr(i). A prefix 
fx of a word lo± ■ ■ ■ CO), is a sub- word of the type io\ ■ ■ ■ uij, with j < k. Let V be the set of 
all prefixes of the words to 1 , . . . ,lo n , together with the N additional symbols ei, . . . , e^, 
which will be regarded as prefixes of the 1-symbol words lo\, . . . , cj± , respectively. We use 
the notation /j, 1 £ V to characterize an arbitrary prefix of lo 1 . Associated to this notation 
there is the function x which extracts the "exponent" of a prefix, meaning that xif 1 ) = * if 
fx is a prefix of uj % . Finally, given a prefix fi of u, we denote by //+ the symbol immediately 
following n in u. This quantity is well-defined, since [i a proper sub-word of lo. 

If we let E^i ^^i x = x + r^i + . . . + r w i , then for £ V we have E^pQi C , and 
hence 

and 

N 

^ = U U 

The same decomposition can be made on each scaled sub-intervals and we thus have a 
recursive tiling of 17. A given point x £ f2 can thus be located if we know for each k 
in which scaled sub-interval it lies. We can show that there exists a sequence of prefixes 
(Mi,/^2> ■ • •) ^ V such that 

oo 
fc=l 

There is a consistency constraint on the sequence, namely 

x(Mfc) = (Mfc+l)+- (21) 

On the possible sequences in "P N , the dynamics of E translates into an odometer-like 
process which can be conveniently represented as a Vershik map, thus justifying the ter- 
minology of "Vershik code". 
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In [18], the periodic sequences of the Vershik code played an important role, firstly 
because it is possible to compute the associated points explicitly, and secondly because all 
the points of the lattice turned out to have periodic or eventually periodic Vershik codes. 
By construction, a point x with periodic Vershik code (/xi . . . pt)°° is an invariant point 
of a contraction 

1 

E^p... E^px = x or x = rp T ^+iP l ( 22 ) 

P 8=0 

where are the translations corresponding to the maps E^, which belong to M. The 
rightmost expression in (|22p may be regarded as an expansion of x in the base p where 
the "digits" r /ii are algebraic numbers (this viewpoint is developed in [19]). Because the 
factors 1 — p T at denominator are not units, in general a periodic point is not an algebraic 
integer, but rather a point of the field K. By the decomposition (Tl2l) . there exists a 
representative £ E S such that x has its orbit on £ + M.. We note that the set £ + M. is 
not invariant under scaling by p, that is, 

pi G £' + M where £' p£ (mod 

Following [18], we notice that, for any positive integer m, the lattice M/m is an O-module 
(see remarks following equation ©), and in particular is invariant under multiplication 
by p. (Alternatively, note that multiplication by p is represented as an integer matrix R 
which preserves the lattices L/m.) With our definition of H, we then see that S D M/m 
has m 3 elements, which limits the number of possible £' after multiplication. Thus, there 
exists an integer t such that p*£ = £ (mod and the smallest such t is called the order 
of £. From equations (|22p . if the point x = (3 + £ with (3 £ A4 and £ 6 H has period T, we 
find that £ — p T ^ G A4, which implies that T must be a multiple of the order of £. 

The period T determines the maximum denominator of the fractional component £ of 
a periodic point, which is given by the determinant of I — R T . Indeed, it is easy to see 
that the characteristic polynomial of I — R T is Pt(1 — t) where pt is the characteristic 
polynomial of R T . Then the eigenvalues of this matrix are of the form 1 — Qf , where Q runs 
through the eigenvalue of R. From the previous section we know that we the eigenvalues 
of R come in reciprocal pairs (£, 1/C), and a straightforward computation shows that 

d r = |det(/-fi r )|=n (Q ^ 

where the product is taken over all the n/2 eigenvalues greater than one. There exists 
then a value of the period above which increases monotonically with T. 

Next we prove a bound for points with eventually periodic code. 

Proposition 9 Let x be a point whose Vershik code is eventually periodic with transient 
t and period T, and let z = <p{x). Then there exists a positive constant C such that 

INI <c||i?||* +T , 

where \ \R\\ is the induced norm of matrices. 
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Proof: On the corresponding lattice, we have 

T-l 



j=0 



and we see that 



T-l 



< ||(7- R^W £ ll#«« + ill < IIC- ^^llmaxj 

3=0 



T-l 

■£W, 

j=0 



and thus 



|z|| < 



1 



ll-RII - 1 

By the discussion above on the eigenvalues of (7 — R T ), we see that 



\\(I-R 



T\-l\ 



max 

k 



1 



max 

k 



1 



max 

k 



s>fc 



i-C 



Now the expression in the max is very small if (k < 1 an d close to one if > 1, it is clear 
that the maximum is attained for the maximum value of |Cfc| ( a t l eas t for sufficiently large 
T), which happens to be ||i?|| and 



||(7 -R T r\ 
We finally have 

\\z\\ < 



max 



a 



( max fc{ii^ji} 



ii^ip 



ii^ir-i 



i + 



ii^ir-i' 



\R\\ - 1 



(||77f -l)<C||i?|p 



Let us now assume that the point z has eventually periodic Vershik code . . . nt(pt+i ■ ■ ■ /J-t+r) 00 ) 
with transient length t and period T. For each positive integer t, we define the mapping 



t-i 



V : L -> L 



z i ^ tp^R-'-tp^Rz = R l z + ^2R 1 v H+1 . 

3=0 



Then V t 1 x has T-periodic Vershik code, and from the above we see that 
INI<q|V t ||||/l|| T <C||ll|| T (||/l||* + 2||lJ||'||t; w+1 ||). 

3=0 

Applying again the geometric sum formula we see that we have a positive constant C' 
such that 

wzWKccwRfmt 

□ 



From this proposition, we know that the number of distinct points z with prefix length 
t and period T which belong to a given layer f + M is of the order of ||i2||( n - 1 )(* +T ). We 
also know that we have up to dr such layers compatible with this specific period T. 
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The number of admissible prefixes of length t is given by the admissibility matrix of 
the graph described by the relation (|2ip . This graph can be described either as a Bratteli 
diagram or as a more conventional connected graph. In any way, it can be shown [21] 
that if a is a primitive substitution, then the admissibility matrix is primitive as well; 
moreover the maximum eigenvalue of this matrix is equal to the maximum eigenvalue 
of M a , the incidence matrix of the substitution. The number of admissible paths in the 
aforementioned graph, as well as all the number of T-cycles, are then of the order of powers 
ofp=\\M a \\. 

We now fix a period T, and consider all the prefix-lengths t. We see that in order that 
all the points fit on all the possible lattices, we must have 

d r ||i2||(»- 1 )( fr +- T ) >C\\M a \\ t+T 

for some constant C. As t increases, this is only possible if ||i?|| n_1 > ||M ff ||. 

We have the following result. 
Proposition 10 If 

WRW"- 1 = \\M a \\, (23) 

then for every period T there is at least one layer M + £ with the property that the set of 
(Vershik) eventually periodic points of period T on the corresponding lattice has positive 
density. 

Proof: The proof follows from the discussion above. Here we mention that the map 
sending each Bratteli code to a point in 0, is not one-to-one; however, it can be shown [21] 
that it is finitely many to one and thus it does not change the result. □ 

If S ^ and n > 2, then the identity (|23|) is never verified. Indeed, in this case, both 
matrices share the same maximum eigenvalue f3, see comments at the end of section 13.11 
This is consistent with the remark we made above, that an orbit must have zero density 
in these cases. 

In the case studied in [18], we had 

ll^ll 2 = ||M ff ||, 

though the result stated here is much less powerful that the results of [18]. 
3.4 Remarks on asymptotic behaviour 

We begin by considering the orbit of the point z which corresponds to the fixed point of 
the substitution. We know that there is an integer j such that the code of z begins with 
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cr k (j) for all k. In this case, we have 



= z + irM^ej 
= z + R k Vj. 

We clearly see that the behaviour here is determined by the largest eigenvalue of R. If 
5 ^ 0, we know that this eigenvalue is (3, and since |cr fc (j)| also scales like f3 k , we have a 
dominant linear behaviour, as expected. 

If S = 0, then the behaviour is dominated by the largest eigenvalue of M a which is 
also part of the spectrum of R. We have 

Tr-D^fcy)! = R k Vj = irM k ej 

where the vector D of discrepancies was defined in ([5]) . The escape rate of the point is given 
by the discrepancy-vector of the point z. There are known results about the discrepancy 
of substitution dynamical systems, including a theorem from Adamczewski [1] which gives 
the growth rate of the maximum component of D^, depending on the properties of the 
spectrum of M a . In particular, we have an upper bound of the form (logk) a k b , with a a 
positive integer and < b < 1. We have this growth when the second largest eigenvalue 
02 of M a is greater than 1 (in modulus). The quantity b has the form 

log |^| 
log/3 

and a + 1 is the multiplicity of 02 

This indicates that a power-law behaviour can only occur when 02 is greater than 1 
in absolute value and has multiplicity one. This eigenvalue then dominates the evolution. 
Indeed if we have 1 1 7ri? i ^.fe (j ) 1 1 1 ~ \@2\ k , since ~ (3 k we are led to 

logjM 

\\nDk\ \ ~ k lo &P . 

In the case studied in [18], direct computation gives 

log I Ail _ 1 
log/3 2 

which leads to positive density orbits on the lattice. 

In general case the computation is less straightforward, but thanks to the recursive 
tiling property we can still derive estimates for the asymptotic behaviour. Let us thus 
consider a point x having Vershik code (m .../*&...). Then 

x g Ef, lP n M n EnpEnf&ixfa) n . . . 

and we can deduce that (see [21] for details) 

x e E^pn^) u{x) = s |Ml| o-(x(^i)) 
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where here s denotes the shift map. To estimate the escape rate of the lattice point 
corresponding to x, we need to know the asymptotic number of occurrences of each letter 
in the words above. Let us call L k the vectors of the number of occurrences of each of the 
letters 1, . . . , N in the k-th word above 

L$ = \s\»M---s M <r(xM)---)\i. 

The number of points in the beginning of the orbit of x thus described is given by the 
length of the word above; since L k is a positive vector, this is given by 

N 

-fell 



Each shift cuts the number of every letter by a finite amount: the number of possible shift 
lengths is equal to the number of possible prefixes, and each one decreases each component 
of L k in a finite numbers of ways. If we call r\ the maximum length of a prefix, then there 
exists not more than N v positive A^-vectors £j of norm not exceeding TV such that 

L 1 = M a e x{fll) - e ai 

I? = M a (M a e x{p2) -£ a2 )- £ ai 

L k = M k e x(pk) -Z-=o MJ ^ 

where 1 < on < N v . As k grows, the escape rate of z under tp is thus given by the 
asymptotic behaviour of 

N k-l 

E L i v * = * Lk = Rkv x(pu) - E Ri ™° j+1 ■ 

i=l j=0 

Since all integer vectors involved have uniformly bounded norm, we have the upper bound 
||^N(*)|| = llvr^H < K ^ ||i?IP < A- "^^ 1 . 

This is not very sharp, yet there is no clear lower bound as one could adjust the £j so as to 
make the L k bounded. The quantity L k is easier to express if x has periodic or eventually 
periodic Vershik code. 

In what follows, we state results only for periodic codes, to which an eventually periodic 
code can be reduced via a finite number of iterations of E. In the previous section, we 
saw that under suitable conditions, these lattice points constitute a "significantly large" 
set, thus justifying our interest in asymptotic estimates. We have the following result. 



Proposition 11 If a point x has periodic Vershik code of period T , then the corresponding 
lattice point z = 4>{x) satisfies the asymptotic estimate 

IT l°gl|fl|l 

U kT {z)\\ ~ fclosll^ll k^OO. 
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Proof: Let x have periodic Vershik code (/ii . . . Ht)°° ■ By formula (|22p . every such point 
is representable as a lattice point. We define 

T-l 
3=0 

for some indices ctj, and 

G(z) = R T z + St. 

The sequence 

nL kT = G k (v xM ) 

is an arithmetic-geometric sequence with respect to k. We use the following lemma 

Lemma 12 If Uk is an arithmetic-geometric sequence, that is 

Uk+i = au k + b 

where a is a square matrix such that I — a is invertible and b is a vector, then 

k 

u k = a k (u Q - I) - I and Uj = (I - ay 1 (I - a k+1 ) (u — l) + kl 

j=0 

with I = (I — a)~ 1 b. 

Proof: We verify that the sequence 

Vk=U k ~l 

is a geometric sequence, indeed 

v k+ i = au k + b- (I - a)~ l b 

= a{u k + a~ l b — a' 1 {I — a)~ l b) 

= a{u k + {fl- 1 {I-a)-ar 1 ){I-a)- 1 b) 

= a{u k - {I - a)~ l b). 

The result then comes by straightforward application of the classic formulae for geometric 
sequences. □ 

We now apply the lemma to our case, with a = R T and b = St- Provided that 
v x{^kT) = (I ~ R T )$T, we obtain 

||V l|LfcT|l (^)|| = ||7rL fcT ||~||i?|| fcT 

Furthermore, if we define 

T-l 

~5 T = Mj>e aj+1 

3=0 
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with the corresponding indices ay and if we call 

G(z) = M T z + St, 

then 

L kT = G fe (e x(/ifcT) ). 
Applying again the argument above with a = Mj and b = 5t yields 

ii^i-iiM.ir 

which leads to 

,rp l°g||fi| 

11^(2)11 ~ fclog||M ff ||. 

□ 

This estimate is very similar to the one obtained for the orbit of zero, a point of constant 
Vershik code. Moreover, as the motion on the lattice is confined to a (n — l)-dimensional 
slab, if 

loell-RII 1 



log P n — 1 

we can expect positive density orbits. 



4 Examples 

In section 13.21 we saw that an algebraic IET of even degree and non-zero drift vector must 
fulfil certain conditions in order to be self-similar: the contraction ratio has to be a root 
of an irreducible self-reciprocal polynomial. A question then arises, naturally, is there any 
such IET? To answer this question, we shall use Rauzy-Veech induction, a powerful tool 
in IET theory, described in the pioneering papers of G. Rauzy and W. Veech [23,25], and 
in many subsequent works. 

Following [25], we let E be an TV-interval IET on the unit interval f2, with irreducible 
permutation ir and lengths vector A. We define 

JV-l N-l 

Io = [0,^2^i) and Ji = [0, K-i®)- 
i=l i=i 

Induction will be on the larger of these two intervals, /. Computing the first-return maps, 
one verifies that there are two cases: 

• |io| < l-^i | an d the induced map has lengths vector Aq(tt)~ 1 A and permutation aovr, 

• |/q| > \h\ and the induced map has lengths vector Ai(n)~ 1 A and permutation a\TT. 
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The permutations a,7r and the matrices Ai are defined as follows: 



irj irj < ttN 

Qokj = \ Tcj + 1 ttN <irj<N 
ttN + 1 irj = N 



CllTTJ 



and 



^o(vr) 




J < TT^N 
j = + 1 

other j 
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V 

where here l n denotes the n-dimensional identity matrix and the first 1 in the last line 
of Aq(tt) occurs at the tt~ 1 N position. Also, if ir~ 1 N = N — 1, the lower-right block of 
Ai(n) contains the single element 1. 



The construction above is naturally associated with a Rauzy graph, whose vertices are 
the irreducible permutations on N symbols. Each vertex has two outgoing edges, labelled 
and 1, ending at the vertices labelled by the permutations of the induced IETs on Iq and 
I\. Each connected component of a Rauzy graph corresponds to a Rauzy class of IETs. 
The Rauzy classes for 3- and 4-interval IETs are shown in figures [U [3j 

An IET on iV intervals is said to satisfy the Keane property [13] if the orbits of its 
discontinuity points are infinite and distinct. Such an IET is both minimal and uniquely 
ergodic; its unique invariant probability measure is Lebesgue measure. If E is an IET on 
iV intervals with irreducible permutation and satisfying the Keane property (this holds 
for almost all length vectors), then, according to a theorem of Rauzy [23], the induced 
map described above is also an IET on iV intervals with irreducible permutation and also 
satisfies the Keane property. 



Rauzy induction provides us with an effective tool for constructing self-similar IETs 
with a prescribed number of intervals. Indeed, a self-similar IET corresponds to a loop 
in the Rauzy graph, with the lengths of the induced map given by a positive eigenvector 
of the iV-dimensional matrix obtained by taking the product of the ^(vr) -1 around the 
loop, and a contraction factor given by the associated eigenvalue. 



4.1 Non-zero drift vector 



For three intervals, the Rauzy diagram has one component and three vertices — see figure 
[TJ From proposition [8] i), the drift vector cannot be zero, for otherwise n would have to 
be 1. But then n must be even, from proposition [8] ii); as we clearly have n < 3, we 
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Figure 1: The Rauzy diagram for 3 intervals 




can only have quadratic cases. This means that the characteristic polynomial p cannot be 
irreducible. 

For four intervals, there are two Rauzy classes, shown in figures [2] and [3j To find self- 
similar IETs, we checked all possible cycles of length at most 16 in a given Rauzy class and 
performed the corresponding products of the matrices Ai(ir). To obtain degree-4 maps, 
we kept only those cycles with primitive product matrices and irreducible characteristic 
polynomials. For such a cycle, initiated at a specific permutation, the lengths of the map 
are the components of the normalized Perron eigenvector of the output matrix. 




Figure 2: The first Rauzy class for 4 intervals. 



Our first example of quartic IET over four intervals comes from the 8-cycle in the graph 
of figureEl determined by the permutations (4213, 4213, 4321, 2431, 3241, 3241, 4321, 4132). 
As in section 13.11 we denote by B the product matrix of this cycle 



B 



( 1 1 1 1 \ 

2 10 

12 2 1 

\ 1 1 1 2 J 



(24) 



with characteristic polynomial 

p(x) = x 4 — 7x 3 + 13x 2 — 7x + 1. 



We see that p(x) is indeed self-reciprocal and we verify that all his roots are real. The 
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Figure 3: The second Rauzy class for 4 intervals. 



smallest eigenvalue of B 

p = i f 7 + V5 - ^38 + Uy/5j , 

corresponds to the normalized eigenvector 

A = (p, 1 - 4p + p 2 , 1 - 4p + hp 2 - p\ -l + 7p- 6p 2 + p 3 ) 

whose entries are positive and sum up to unity. These data, together with the permutation 
(4213) define an IET E, with module M = Z[p], for which we choose the canonical basis 
Vi = p 4-1 , i = 1, ... ,4. One verifies that the induced map on the first interval Oi is just 
the original map E rescaled by a factor p. Figure 0] shows the action of the map. 
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Figure 4: The interval exchange transformation E 



The following matrix transforms the lengths into the translations vector r = A • A: 



A = 
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1 \ 




-1 





1 







-1 


-1 










-1 








/ 



The IET E is conjugated to its first return in the interval f^i via the substitution 

'I — ► 143 
I 2 — ► 143223 

3 — ► 14323 

4 — > 1443. 

We note that the transpose of the inverse of B in (|24p is the incidence matrix of a. This 
identification is generally true for Rauzy induction, due to the identification of the zth 
substitution string with the first-return itinerary of the ith IET sub-interval. 
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Finally, we verify that the drift vector is non-zero. The integer translation vectors 
Vi = 4>(Ti) — see equation ([6]) — are 

(1,-1,0,0), (1,-5,5,-1), (-1,3,-1,0), (0,-1,0,0), 

so that 

4 

S = Ai«j = (p- 4p 2 - p 3 , -1 + 16p 2 - V, 4 - 16/9 + p 3 , -1 + Ap - p 2 ) 

i=i 

is non-zero (since p is an algebraic number of degree four). 

The cycle of length 8 containing E is the only one with a quartic irreducible polynomial. 
The other members of the cycle give 7 other IETs which renormalize with a scale factor 
p, but with different permutations and sets of widths in the ring 7L\p\. The example E 
has the nice property that the rescaled interval is just fij, a feature not shared by all of 
its cycle-mates. For example, the cycle contains another IET with permutation (4213) for 
which the induced map again acts on the interval (0,p),p ~ 0.227777, but the latter is 
properly contained in Oi = [0, p), p = —1 + 7p — 6p 2 + p 3 0.294963. 

Interestingly, for longer cycles we also have relatively few irreducible quartic polyno- 
mials (hence rings). The table below gives a summary of the numbers of suitable Rauzy 
cycles and associated polynomials for cycle length up to 16. 



cycle length 


U Rauzy cycles 


tt polys 
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80 


27 
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148 


50 


14 


360 


108 


15 


798 


211 


16 


1657 


452 



4.2 Zero drift-vector 

Here we are interested in IET belonging to the Arnoux-Rauzy family [6] to which belongs 
the Arnoux-Yoccoz example [4,7] studied in great detail in [18,19]. Originally, this family 
was found as a way to construct a sequence of complexity 2n + 1. The construction is the 
following (we keep the original notation from [6]). 

Let a, (3, 7 be three positive real numbers such that a > (3 + 7. We define a six- 
interval-exchange map / over the circle of length a + (3 + 7 as follows. First, we divide 
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in two halves the three intervals I a , Is, J 7 with respective lengths a, (3, 7 and we swap 
the halves preserving the orientation. Then we rotate the whole circle by it. It can then 
be shown than the induced map on f(I a ) has the same structure over three intervals of 
respective lengths a' = a — [5 — 7, = /3, 7' = 7. We can iterate this process: at each 
step we induce on the largest interval, and at each step one of the intervals of the induced 
map is larger than the sum of the other two. We can see that this process has the same 
flavour as Rauzy induction. If we define the three matrices 

/ 1 1 1 \ / 1 \ / 1 \ 

Mi = 1 M 2 = 1 1 1 M 3 = 1 

\ 1 / \ 1 / \ 1 1 1 / 

then it can be shown that for any k the product M1M2M3 occurs infinitely many times 
in the possible sequences of induction. That is, for each k, the Perron eigenvector of 
each such product will give the lengths of a renormalizable map E k over the circle. By 
construction, all members of this family have zero Sah-Arnoux-Fathi invariant, and hence 
zero drift vector (proposition [5]). 

The case k = 1 is the Arnoux-Yoccoz case with three intervals, with characteristic poly- 
nomial p(x) = x 3 — 7x 2 + 5x — 1. Its unique real root is the scaling constant p = A~ 3 , where 
A is the so-called tribonacci number [22, p. 233]. The lengths are (a,/3, 7) = (A~ 2 ,A _1 ,1). 
This example has very interesting properties, which include finite decomposition. 
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Figure 5: The interval exchange Ei- 



For k > 1, the maps maps E k are quite similar to the Arnoux-Yoccoz map. We will 
present some computations for the cases k = 2,3,4,5,6 (see figure for k = 2), which 
show that they are likely to share the finite decomposition property. 

The module of E^ is the half-integer lattices Ai = Z[Afc]/2, where A& is the smallest 
real roots of the polynomial given by 

f k {x) = x 3 - (k + 4)x 2 + (3k + 4)x - 1, k > 1. 

One verifies that for all A; > 1, the polynomial ft has no integer roots and hence is 
irreducible. The (non-normalized) vector of lengths of the partition is 

A = (A fe - Xl/2, A fc - A^/2, 1/2 - 3A fc /2 + \\/2, 1/2 - 3A fc /2 + A^/2, 1/2, A fc /2, 1/2 - A fc /2), 

and the corresponding translation vectors are 

+ Afc — A 2 , 2 — 3Afc + A^, 3 — A\ k + A 2 , 1 + 2A^ — X 2 ., A& — 1, 1 — A^, A& — 3). 

One can then verify directly that the drift vector is indeed zero. The maps E k are not 
renormalizable but in each cases — as in [18] — the induced map in the first interval leads 
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to a renormalizable IET on seven intervals. We shall carry out the computations directly 
with the maps E^. 

In order to investigate the finite decomposition property, we will perform a fairly 
straightforward algorithm implemented in C. For each k we will first select the slab of the 
3-dimensional half-integer lattice consisting of the points belonging to the interval [0, 2— A*.) 
which lie inside a cube of prescribed size. Then we iterate forward and backward all the 
points in a small cube around the origin to check how many points in the big cube slab are 
attained. This method is fast, easy to implement, and applicable to any algebraic IET. 
It is however memory-intensive, which constrained the maximum linear size of the cube 
we were able to consider to about 600 points. For the sake of speed, when mapping a 
lattice point into the intervals we use standard double precision (as opposed to arbitrary 
precision). It should be noted however that there is no error propagation in the dynamics, 
as each point is represented in vector form using the Z-module basis. 

To represent lattice points, we chose the basis (1, A&, \\) for Z[A]. It should be noted 
that this choice is not necessarily the best in all cases. Indeed, orbits are often highly 
"anisotropic" , escaping much faster in a specific direction than in others. As a consequence, 
covering a significant part of the cube around the origin could be time-consuming. While 
more favourable bases have been found in special cases, no general optimal family has 
been found. For instance, for k = 1 we have the same case as in [18] and yet the orbits are 
not distributed as uniformly as in this paper. In fact, Ai is the third-root of the tribonacci 
number used as a basis in [18], switching to this basis gives more uniform orbits. This 
process does not generalize to greater k though. 

Nevertheless, on all the cases we have investigated we see that a very small set of initial 
conditions is required to cover at least a significant part of the lattice around zero. We 
define 

Cd = {\(p,q,r) e hi? : ±(p + q \ k + r \l) G [0,2 -A fc ), max{|p|, |g|, |r|} < 2d}. 

In table[H we summarize the results of the investigations. For each k, we give the maximum 
size D of Cd that we were able to fill with backward and forward iterates of points in an 
initial region Cd of size d. We also give an upper bound on the number T of iterations 
required. 
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D 


d 


T 


2 


400 


1 


3 x 10 c 


3 


400 


3 


6 x 10 7 


4 


400 


1 


2.6 x 10 8 


5 


400 


1 


9 x 10 8 


6 


200 


3 


1 x 10 9 



Table 1: Data for lattice- filling experiments. 

For each k, the substitutions au are explicitly computable and we can use propositions 
[T0l and [TT] to give us additional evidence regarding the finite-decomposition hypothesis. 
Indeed the latter proposition gives us asymptotic estimates on the escape rate of orbits 
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of the lattices having periodic Vershik code, while the former gives a condition for these 
orbits to be fairly common. Both proposition rely on the value of the following quantity 



_ log||i?|| 
V ~ logHM^H 

which gives the exponent of the escape rate of the orbit of a point having periodic Vershik 
code. These points have density one in at least one coset £ + «M if v = l/(n — 1). In the 
following table we give approximations of the values of v for several values of k 



k 


1 


2 


3 


4 


5 


6 


7 


V 


0.5 


0.5 


0.5 


0.546715 


0.595958 


0.623202 


0.642502 



We see that for k = 1,2, 3, we can deduce from the proposition mentioned above that 
the points having periodic Vershik code (a set which have positive density) have positive 
density, which is another indication that E 2 , E% might have the finite decomposition prop- 
erty. The data for the other values of k are inconclusive though. The following table shows 
values of v for increasing values of k. The value of v appears to increases slowly toward 
unity. 



k 


17 


30 


50 


100 


150 


200 


500 


V 


0.721655 


0.756235 


0.780871 


0.8074 


0.820165 


0.828247 


0.849766 



This can be explained if we look at the partitions of the maps E^ for increasing k. 
As displayed in figure ([6]), we see that the intervals carrying the "rotation- flip" dynamics 
tends to fill the whole interval, leading to a global dynamics approaching a simple non-zero 
drift dynamics, for which v = 1. 



4.3 Self-similar E 2 model 



We conclude our discussion of zero-drift-vector examples with a more detailed look at a 
self-similar version of the map E 2 introduced in section 14.21 The particular IET, which 
we shall call E 2 , is a member of a 29- vertex cycle in the same 294- vertex Rauzy class 
of 7-interval IETs as the self-similar Arnouz-Yoccoz model of [18]. Its permutation is 
(5462731), and the 7-dimensional cyclic product matrix is 
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equation 


(ED 


))) 









P(x) = O 3 - 6x 2 + lOx - l)(x 3 - 10x 2 + 6x - 1)(1 - x). 
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k = 2 



k = 3 



k=lO 



k = 20 



Figure 6: The interval exchanges E^ for increasing k, showing that three intervals tends 
to vanish, and the dynamics approaches a rotation with flips. 

The scaling constant p ~ 0.106711 is the real root of the polynomial 

p(x) = x 3 — 6x 2 + Wx — 1. 
One verifies that = p~ 1 , the real root of p(x) = x 3 — Wx 2 + 6x — 1, is a Pisot number. 

The 29 members of the Rauzy cycle containing E\ are an self-similar IETs with the 
same scale factor p. The map E\ was selected for special consideration because of several 
additional simplifying features. In particular, its translation module is A4 = 1i[p], its 
rightmost interval IV is of length p, and the induced IET on Q7 is just the original map 
rescaled by p. 

The vector A of interval lengths is the positive 7-vector associated with p in the diag- 
onalization of B^ 1 , normalized to total length unity. Specifically, choosing the canonical 
basis (l,p, p 2 ) for M = Z[p], we find 

A = (1 — 5/9 + 2p 2 , -1 + 10/9 - 3p 2 , 1 - 9p + 3/5 2 , p - p 2 , -1 + lip - Ap 2 , l-9p + 3p 2 , p) 

with corresponding integers translations Vi = 4>{ji), i = I, . . . ,7, 

( Vi ) = ((0, 3, -1), (0, -2, 0), (1, -7, 2), (-1, 5, -2), (1, -8, 3), (0, -6, 2), (-1, 1,0)) . 

The IET is illustrated in figure [71 Calculating ^ AjUj, we verify directly that the drift 
vector vanishes. 

Iterating the map E%, we can examine the first-return orbits of the 7 intervals of the 
induced map on ^7 and establish the recursive tiling property, with associated substitution 
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2 3 4 5 6 7 



7 4 6 2 



3 5 



Figure 7: The self-similar E 2 map. 



rule 



1 — ► 7114115 

2 — ► 7114121361361214115 

3 — ► 711412136136135 
rr : { 4 — ► 71141214115 

5 — ► 71156136135 

6 — ► 711561361361214115 

7 — ► 7115. 



We obtain the following result: 



Proposition 13 We consider the mapping E^. The following holds 



i) Let x be a point of the unit interval and let fJ,(x) be its Vershik code. Then x £ 
Q(p) n [0, 1) if and only if p,(x) is eventually periodic. 

ii) The restriction of the module M to the unit interval consists of the (forward and 
backward) orbits of the discontinuity points. 

Hi) The map has the finite decomposition property. Specifically, there is a positive integer 
M such that for each £ £ K the restriction of the layer £ + M. to the unit interval 
decomposes into at most M orbits. 



We shall not provide a detailed proof of the proposition above, since it is a straight- 
forward adaptation to the present model of the proof provided in section 3 of [18] for the 
Arnoux-Yoccoz cubic model. Let us merely summarize the main ideas. We represent a 
point x in Q(p) D [0, 1), as x = £ + z, with ( £ S — see equation (fl~3j) — and z G 
Specifically, using the basis (1, p, p 2 ) for Q(p) we have 

£ = r + r lP 2 + r 2 p 2 n £ Q n [0, 1) 
z = mo + m\p + m,2p 2 mi G Z. 

The restriction of x to the unit interval leads to the condition 

m = - + m\p + m 2 p 2 \ ■ 

Any point on the forward or backward orbit of x can be represented in the same way, 
with the same £. Since each point is completely determined by an integer pair (1711,1712), 
it is natural to associate the orbit with a walk on the lattice = Z 2 (with reference to 
equation 0, here we have n = 3 and 6=1). 
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To prove the finite decomposition property of proposition [131 we must show that the 
number of lattice orbits in is bounded uniformly in £. The strategy is to show that 
every lattice orbit visits a finite ^-independent core region containing no more than a finite 
number M of lattice points. The existence of such a region is also the key to proving the 
second statement in proposition I13t one needs only to verify that the lattice orbits of the 
discontinuity points visit all of the points of the core region of Lg. 

To complete the proof, one shows that for each £ 6 S, some kth iterate of the left-shift 
map for the Vershik code is represented in by a map j k with the following property: 
every (mi, 1112) £ has a forward 7^ orbit which is eventually trapped in a finite, £- 
independent core region. This orbit must eventually settle down in a limit cycle of some 
finite period t. Thus the original point x has an eventually periodic Vershik code of period 
kt, thereby establishing the "only if part of proposition [13] i); while the "if part follows 
at once from (|22p . Since the periodic tail of this code coincides with that of a point in the 
limit cycle in the core region, the two points lie on the same E% orbit, and we have the 
remaining element in the proof of the finite decomposition property as well. 

It is important to point out that the proof above relies essentially on the self-similarity 
of the IET, and on the Pisot property of the reciprocal of the scale factor. In particular, it 
is the latter which ensures a contraction to the core region. In [18], additional properties of 
the specific IET were used to characterize in some detail the lattice orbits in the various L^. 
The orbits of the self-similar Arnoux-Yoccoz model were found to be non-crossing, filling 
out with unit density disjoint sectors of the lattices. It turns out that the lattice orbits of 
the E% model do not have these features. This is already true for the discontinuity-point 
orbits, for which the lattice is partitioned into three sectors, each filled out by a pair of 
mutually crossing orbits, neither of which is locally lattice-filling. The situation near the 
origin is shown in figure [8] 
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